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Abstract. I revisit a longstanding question in quantum optics; When is the rotating wave 
approximation justified? In terms of the Jaynes-Cummings and Rabi models I demonstrate that 
the approximation in general breaks down in the adiabatic limit regardless of system parameters. 
This is explicitly shown by comparing Berry phases of the two models, where it is found that this 
geometrical phase is strictly zero in the Rabi model contrary to the non-trivial Berry phase of the 
Jaynes-Cummings model. The source of this surprising result is traced back to different topologies in 
the two models. 



1. Introduction 

The rotating wave approximation (RWA) is taught in any graduate course on quantum optics [JJ. One 
of the more familiar examples were the RWA has proven very efficient is for the Jaynes-Cummings 
(JC) model describing the interaction between a single two-level atom and one quantized resonator 
mode [U [2] . For moderate field amplitudes and quasi-resonant interaction between the field and atomic 
transition frequencies, the approximation is taken to be justified whenever the atom-field coupling 
is small in comparison to the mode frequency. This has been verified by numerous cavity QED 
experiments [3], where the JC model correctly predicts the measured quantities. 

It is not until very recently that the applicability of the RWA in the adiabatic limit was 
questioned |3J. We know that for very long evolution times, small perturbations can accumulate 
and greatly affect the results. Of particular interest is when the perturbation alters the symmetry 
of Hamiltonian, or even its topology. In this work I show how this is exactly the case for the JC model. 
At a mean-field level, the adiabatic energy surfaces in phase space change topology when the RWA is 
applied. In particular, the surfaces of the JC model possesses a conical intersection (CI) absent for the 
Rabi model (JC model including the counter rotating terms). Such a CI is the origin of non-trivial Berry 
phases [5] , and as a result one concludes that the commonly considered Berry phase must strictly vanish 
in the Rabi model while it is non-trivial for the JC model. Employing numerical diagonalization, going 
beyond the above mean-field arguments, it is also demonstrated that the so called "vacuum induced 
Berry phase" vanishes. 

2. Adiabatic concepts 

Adiabaticity has many appearances in quantum mechanics. The original idea dating back to 1928 [B] 
and given in most advanced text books on quantum mechanics considers an explicitly time-dependent 
Hamiltonian H(t), and the adiabatic theorem states that a system prepared in an instantaneous 
eigenstate of H(t) remains in this state provided the change of H(t) is "slow enough" and there 
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is an energy gap to nearby eigenenergies. While the adiabatic theorem deals with strictly time- 
dependent systems, the notion of adiabaticity can be efficiently translated to other situations where 
the time-dependence is replaced by an intrinsic evolution. The system evolution can, in many cases, be 
characterized by various time-scales. Whenever these scales are very different, it become legitimate to 
perform an adiabatic separation or elimination of the fast degrees-of-freedom. This is the basis of the 
Born-Oppenheimer approximation (HE]- 

Consider a system with two well separated characteristic time-scales. Letting R and r denote 
dynamical variables for the slow and fast degrees-of-freedom respectively, we write the Hamiltonian 
H(R,r). Within the BOA we assume the fast degrees-of-freedom to instantaneously follow the motion 
of the slow degrees-of-freedom. Thus, freezing R we can solve the eigenvalue problem H (R, r)\Q(R, r)) = 
V a d(R)\&(R,r)), where |$(i?,r)) is the adiabatic eigenstate and V a d{R) the adiabatic potentials. The 
adiabatic eigenstates (labeled them with some quantum number(s) n) form a complete basis for the full 
system state |^(-R, r)). The BOA consists in neglecting the (" non-adiabatic" ) couplings A(R, r) between 
these adiabatic eigenstates, and the corresponding coefficients 4> n {R) in an expansion of \^>(R,r)) are 
given by the eigenvalue problem H(R)<f)(R) = T + V a d(R) 4> n (R) — E(j) n (R), where we symbolize the 

"kinetic" energy operator for the slow degrees-of-freedom with T. 

For this work, where we study a spin 1 /2 particle coupled to a single boson mode, the slow degree- 
of-frccdom is the boson mode representing an electromagnetic resonator mode and the internal two-level 
structure of the particle is the fast degree-of-freedom. Typically, the characteristic instantaneous time- 
scales depend on the full state \^(R,r)) and if there is no clear separation in the time-scales the BOA 
breaks down, which become evident in the fact that the non-adiabatic couplings A(R,r) are no longer 
negligible. Note, however, that the absolute size of \A(R, r)\ does not serve as a measure of adiabaticity, 
but it must be put in relation to the actual dynamics [8]. A(R, r) is the so-called Mead-Berry connection 
and has the properties of a gauge potential and is in general on matrix form [5]. The gauge freedom 
lies in the ambiguity of an overall (independent) phase of the adiabatic states \&(R, r)). 



3. Jaynes-Cummings vs. Rabi 

3.1. Breakdown of the RWA in the ultrastrong coupling regime 

For the last three decades, the JC model has served as a workhorse in cavity QED 2J, explaining and 
predicting many of the cavity QED experiments to date j3] . With the development of circuit QED [9] , 
the validity of the RWA applied in the derivation of the JC model has been thoroughly explored, see 
for example [TO]. The JC model is derived from the Rabi Hamiltonian (% — 1) 

Hr = loo) a + + ffV^ (a) + a) a x (1) 

by neglecting the counter rotating terms a* a + and aa~ corresponding to virtual light-matter excitations, 
i.e. 

Hjc = uja^a + —a z + gy/2 (a) a~ + <7 + a) . (2) 

Here uj and f2 are the respective resonator and qubit (atom) transition frequencies, g the effective 
qubit-light coupling, (a) is the photon creation (annihilation) operator, and the u-operators are the 
standard Pauli matrices acting on the internal qubit states |1) and |2). In this setting, the justification 
of the RWA is typically given by u > g [11]. While cavity QED experiments operate in the regime 
10~ 7 < q/ui < 10 -5 [3], the first generation circuit QED realizations reached g/ui ~ 10~ 2 [12 , and 
more recently g/ui ~ 0.1 has been achieved and as a consequence the Bloch-Siegert shift measured [T3]. 
For g ^> k, r where k is the characteristic resonator decay rate and T the decay of the qubit, the short 
time evolution is predominantly coherent and the system is said to work in the strong coupling regime. 
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The works of ref. |13j operate instead in the ultrastrong coupling regime in which g/u) > 0.1. The eleven 
lowest eigenenergies for the Rabi and JC Hamiltonians as a function of the coupling strength are shown 
in figure [I] Throughout this work we use lo = 1, i.e. frequencies (energies) are scaled in terms of the 
photon energy. In agreement with the above mentioned experiment, at around g ~ 0.1 the Bloch-Siegert 
shift becomes clearly visible with a relative error ~ 0.1 %. Further, at g = y/2 the JC spectrum shows 
a crossing and the eigenstate |0, 1) (vacuum and deexcited qubit) is no longer the ground state. For 
the same coupling, when including N identical qubits, hence considering the Dicke model |14) . the Rabi 
model possesses a thermal |15j as well as quantum phase transition 16J from a normal to a superradiant 
phase. 




Figure 1. The eleven first eigenvalues E v (g) (y = 1, 2,..., 11) for the Rabi (black) and JC (green) 
models. The frequencies are ui = Q = 1. 

The conservation of excitations N — n + \a z (h = a'a) of the JC model results in a continuous 
U(l) symmetry and thereby integrability. Together with the "empty" eigenstate |0, 1), the eigenstates 
are the dressed states, 

|$+) = cos(0)|n-l,2>+ S in(0)|ri,l}, 

(3) 

|$-) = cos(0)|n- 1,2) +sin(0)|ra,l), 
where tan(20) = 2gy/2n/A with n a positive integer and A = Q — u the qubit-field detuning. The 
counter rotating terms break the U(l) symmetry and the Rabi model has instead a discrete Z 2 parity 
symmetry. It was not until recently that it was demonstrated that also the Rabi model is (quantum) 
integrable [17] . The analytic solutions, however, are not on simple closed form; both the spectrum and 
the eigenstates are given as sums of complicated functions. Thus, the search for simple approximate 
expressions for the solutions of the Rabi model has continued, see among others [18] . 



3.2. Mean- field Born-Oppenheimer approach 

Let us now employ the adiabatic ideas of the previous section to the JC and Rabi models. We have 
the degree-of-freedom for the boson mode and that of the qubit. As already pointed out, the validity 
of the BOA will depend on the actual state we consider. We begin by rewriting the Hamiltonians in a 
quadrature representation 

Hjc = ^&z+g (xa x + pay) (4) 

and 

Hr = oj f — + y j + —b z + 2gxa x . (5) 
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Here we give the JC Hamiltonian in the interaction picture, and the quadratures are 

1 i 
x= -=(a + ff) , p=-j=(a-a < ). (6) 

In the terminology of the previous section, R represents x and p while r serves as the a a (a — x, y, z). 
Applying the BOA, we diagonalize (fast evolving) spin degrees-of-freedom under the assumption of 
commutability of x and p. Within this approximation we identify the resulting "diagonal" Hamiltonian 
as semiclassical energy surfaces E±(x,p), explicitly given by 




Figure 2. The semiclassical energy surfaces for the JC (a) and Rabi (b) models. The JC model shows 
a CI absent in the Rabi model. 

The semiclassical energy surfaces of figure [2] make clear that the presence or absence of counter 
rotating terms change the general structure. For the JC model there is a CI degeneracy at the origin 
for vanishing detuning A = 0. In the Rabi model there is not a single point degeneracy but a seam 
along the p-axis. The energy gap SE = E + (x,p) — E_(x,p) sets the characteristic time-scale for the 
spin, and it follows that in general large photon field state implies justification of the BOA, i.e. for 
n = p 2 /2 + x 2 /2 ^> 0. The exception occurs for x » in the Rabi model where the gap is simply 
SE = ft. 

4. Adiabaticity and Berry phases in cavity/circuit QED 

Symmetry, topology and geometry play important roles in physics. For example, two Hamiltonians can 
have nearby identical spectrum, but different symmetries or topologies which may have far-reaching 
consequences. In recent years these aspects have seen especially renewed interests, for example in 
topological quantum computing |19j and topological insulators [20] - Other known examples worth 
mentioning and related to this work include: 

(i) Geometric Berry phases. When some state vector ^(i?)), being an instantaneous non-degenerate 
gaped eigenstate of the Hamiltonian, is taken adiabatically around a closed loop C in parameter 
space R, it acquires an overall additional phase 7c to the dynamical one |21) . Hence, the parameters 
R are indirectly time-dependent. This phase has a geometric origin and depends on the loop C, 
but not on the actual time it takes to enclose the loop. Explicitly it reads [3T] 

7c = / dR- A(R) = %l dR- (V(R)\V R \y(R)), (8) 
Jc Jc 
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where we have reintroduced the Mead-Berry connection A(R) [5]. As for the Mead-Berry connection 
for the Born-Oppenheimer treatment, A(R) is again comprising the non-adiabatic couplings. The 
Berry phase Q is defined in the strict adiabatic limit, where we must have that the non-adiabatic 
couplings vanish. Nevertheless, even though the non-adiabatic couplings are zero we may have 
7c ^ (or any multiple of 2tt). This reflects the geometric property of the Berry phase; 7c must 
depend on the geometry of the energy surfaces in parameter space. 

(ii) Reduction of symmetry. We already touched upon the fact that the continuous U(l) symmetry of 
the JC model is reduced to a discrete Zi symmetry when the counter rotating terms are included. 
This implies that the corresponding classical models (obtained in a mean-field coherent state ansatz) 
are either integrable (JC) or non-integrable (Rabi). As a consequence, the JC model does not show 
classical chaos while the Rabi one does [22] , 

(iii) Non-commutability between semiclassical and adiabatic limits. The classical limit of quantum 
mechanics is often taken as H — > 0. The adiabatic limit is assumed to be that for which any 
characteristic time-scale T goes to infinity, v = T _1 — > 0. It has been argued that these two limits 
are equivalent |23j . This was later proved wrong by considering a counter example in terms of a 
Landau-Zener problem [24 . More precisely, the order in letting Ti and v approach zero may render 
qualitatively different results. In a later reference [25 , where a non-linear Landau-Zener model was 
analyzed, it was demonstrated that the non-cummatibility between the two limits appeared when 
the instantaneous eigenergies change topology. 

(iv) Secular approximation. Utilizing several approximations is the standard procedure when deriving 
master equation describing the reduced dynamics of a system weakly coupled to an environment. 
The Born approximation neglects correlations between the system and the bath, the Markovian 
approximation disregards memory-effects of the bath, and finally the secular approximation is a 
version of the RWA |26j . The application of the secular approximation casts the master equation 
into a Lindblad form known to preserve many physical properties |27j . Despite this appealing form, 
it has been shown that the secular approximation can result in unphysical results [28] . 

Having argued that the topology of the energy surfaces in parameter space can affect system 
properties greatly, we return to the JC and Rabi models. Already figure [2] signals a qualitative difference 
between the two models. The semiclassical energy surfaces of the JC model possesses a CI. It has long 
been known, both in molecular/chemical physics [5] [29] and in condensed matter theories 130,, that CI's 
render the system dynamics, and it is especially recognized that encircling a CI adiabatically results in 
a non-trivial geometric Berry phase. At the same time, it is known that without CI's there are only 
trivial Berry phases [5] ■ It is thereby expected that the JC and Rabi models give different Berry phases 
when encircling the origin (x,p) — (0,0). 

Our two models are time-independent and to compass the origin we apply the unitary operator 
U(4>) = exp [—incj)] to our Hamiltonians [31 j . and think of </> as an external controllable parameter that 
is varied adiabatically 4> : — > 2tt. Experimentally, such a change can be accomplished via an external 
driving field with adjustable phase. The effect of U(4>) is to rotate the state in phase space, i.e. bringing 
about the encircling of the origin. 

Employing the BOA we can directly calculate the resulting Berry phase for the JC model [32] 



where \$±(4>)) are the two adiabatic eigenstates of the JC model, and p is related to the field amplitude 
as p 2 /2 = p 2 /2 + x 2 /2 = n + 1/2. The above result relies an the BOA where we neglect corrections 
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from the non-commutability between x and p. If we compare to the exact result using the eigenstates 
Q one finds an 1 /2 error [3T] : 

i [ i c) =±Ai- , . A/2 i. 



In the large field limit, where the BOA becomes exact, the two results agree. We may rephrase the 
mean-field Born-Oppenheimer analysis by directly study the unitary transformed Hamiltonian 

H'j (4>) = U(<P)H JC U- 1 (<P) 

(11) 

= + [(cos 4> x + sin (ftp) cj x + (sin i — cos <^p) <7j,] . 

Introducing an effective magnetic field B(0) = (B x , B y , B z ) — (cos a; + sin </3p, sin 0i — cos <frp, A/2) 
we can write H'j C (<j)) = B(0) • <r mimicking a spin 1/2 particle in a magnetic field. Since all three field 
components B x , B y and B z are in general non-zero there is no way the adiabatic eigenstates |$±(y>)) 
can become strictly real-valued and jjc will typically be non-trivial. 

Turning now to the Rabi model we sec that the transformed Hamiltonian reads 

( p 1 x 2 \ 

H' R = u f — + — + — a z + 2g (cos^x - sin<pp)a x . (12) 

The effective y-component magnetic field B y — implying that there is a gauge freedom to choose the 
adiabatic eigenstates real-valued with a vanishing Berry phase as a result. As we already pointed out, 
this result follows also directly from the semiclassical energy surfaces of figure [2] since the Rabi surfaces 
E { ±\x,p) lack any CI's 0. 
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Figure 3. Numerically calculated Berry phase for the first excited state of the Rabi model. The lines 
correspond to different coupling strengths, g = 0.001 (green), g = 0.01 (black), g = 0.1 (red) and 
g = 1 (blue). Note that independent of g, the Berry phase always vanishes after closing the loop C at 
4> = 2tt. 

The absence of a Berry phase in the Rabi model as declared above relies on mean-field 
approximations, and it is not evident that these findings hold true also for small field amplitudes. 
Nevertheless, the change in topology in the energy surfaces between the two models suggests that also 
in the weak field limit the Berry phase should disappear. We verify this statement by numerically 
calculating the Berry phase for the Rabi model. The eigenstates of the model are obtained by 
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diagonalizing the Hamiltonian H' R ((f>) expressed in the bare basis {|n, 1), |n, 2}} where n is a non- 
negative integer truncated at some large value (~ 500). The Berry phase ([8]), ^ is then calculated 
by numerically evaluating the integral for <fi : — >■ 2tt. We consider the first excited state since the 
groundstate has a vanishing Berry phase also for the JC model [31]. In the JC model, this state can 
be thought of as initially the atom is excited while the cavity is in vacuum. The results are displayed 
in figure [3] showing how the Berry phase varies as <\> : — > 2tt for different couplings g. What is clear is 
that regardless of g, the phase is strictly zero when closing the loop (</> = 2tt). This demonstrates how 
"vacuum induced Berry phases" do not exist when including the counter rotating terms. We may note 
that g — 0.001 is in between cavity QED and typical circuit QED experiments, g — 0.01 is in agreement 
with the circuit QED experiments [T2], while g = 0.1 are in the ultrastrong coupling regime recently 
demonstrated in [13j . 

5. Concluding remarks 

The results presented in this work can be taken as being in conflict with earlier studies |31U33j . However, 
all of them start from a rotating wave Hamiltonian and as demonstrated here, it is then not surprising 
that the Berry phase is non-trivial. Interestingly, long before the work by Fuentes-Guridi et al. the 
vacuum Berry phase of the JC model had been considered in reference }34| with contradicting results, 
i.e. they claim that this phase vanishes also in the JC model. It should be noted though, that the 
employed WKB method of [M] relies on semiclassical arguments. Despite the fact that the Berry phase 
seems to disappear in the Rabi model, I am not saying that irrelevant Berry phases are not achievable in 
cavity/circuit QED systems. Other system configurations can indeed generate topologically non-trivial 
energy surfaces, for example in bimodel set-ups [35] . 

Another relevant aspect concerns the accuracy of the JC or Rabi models to describe certain physical 
systems. Already performing the single mode and two-level approximations utilize sort of adiabatic 
elimination schemes, and it is not clear whether such terms can as well influence the topology and 
thereby alter results in the adiabatic limit. It is generally believed that the Rabi model presents a more 
precise description of current cavity and circuit QED experiments than that of the JC model. This 
seems as a trivial statement, but the issue is indeed more subtle than believed [33]. 

Maybe the most striking outcome of this work relates to the experiments that claim to have 
measured the Berry phase in circuit QED [37]. The JC results agree well with these experimentally 
extracted phases. One should remember, though, that the Berry phase is a strict adiabatic result, and 
no experiment can ever reach this limit, mainly due to definite coherence times but also possibly due 
to finite lifetimes of the experimentalists. In this respect, to be mathematically correct their results are 
rather measurements of the geometric Aharonov- Anandan phase [35] , for which my results do not apply. 
The natural question then follows; What are the time-scales for the effects of the counter rotating terms 
to become visible? Numerically answer this question would be very demanding. Almost all numerical 
methods rely on some sort of truncation, and as the characteristic time-scale grows very large it will 
most likely imply that the truncation limit has to be pushed further and further in order to keep the 
numerical accuracy, implying increasing computational power. Related to this issue, we note that there 
are some recent works contradicting my results |39j . Both these use series expansions for approximating 
the solutions of the Rabi model, and I argue throughout this work that in the adiabatic limit one must 
be very careful with any neglected small terms, and any truncation of the series expansion must thereby 
be carefully justified. 
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